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The "overlapping-cycles shuffle" mixes a deck of n cards by mov- 
ing either the nth card or the (n — fc)th card to the top of the deck, 
with probability half each. We determine the spectral gap for the 
location of a single card, which, as a function of k and n, has sur- 
prising behavior. For example, suppose k is the closest integer to an 
for a fixed real a £ (0,1). Then for rational a the spectral gap is 
0(n -2 ), while for poorly approximable irrational numbers a, such as 
the reciprocal of the golden ratio, the spectral gap is 0(n~ 3/ ' 2 ). 

1. Introduction. The overlapping-cycles shuffle has two positive integer 
parameters n and k £ [l,n), and mixes a deck of n cards by flipping a fair 
coin to decide whether to move the nth card or the (n — k)th card to the 
top of the deck. The case where k = 1 was first studied by Hildebrand [6], 
who showed that the shuffle mixes within 0(n 3 log n) time, and Wilson [8] 
who showed that it takes at least (1 — o(l))/(87r 2 )n 3 logn shuffles to mix. 
Jonasson [7] first studied the overlapping-cycles shuffle that we consider 
here, obtaining a lower bound of order ra 3 //c 2 logn, but for this shuffle there 
was no matching upper bound, and indeed there cannot be one, since, as 
Jonasson showed, when k = n/2 the shuffle takes at least order n 2 time to 
mix. 

We would like to understand the mixing time of this overlapping-cycle 
shuffle, but the behavior is sufficiently complicated that we settle for the 
more modest goal of understanding how a single card behaves under this 
shuffle, and in particular, what the spectral gap is for the Markov chain 
induced by a single card. The behavior of the spectral gap is a surprisingly 
complicated function of n and k (see Figures 1 and 2), which depends to a 
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Fig. 1. The relaxation time as a function of n and k, for 1 < k < n < 200. 

great extent on the Diophantine approximations of k/n. For example, sup- 
pose a is a fixed number between and 1, and k = [an] , that is, the closest 
integer to an. Then the inverse spectral gap (also known as the relaxation 
time) is of order n 2 when a is rational, but when a is a badly approximable 
irrational number, such as the golden ratio minus 1, the relaxation time is 
of order ?i 3 / 2 . 

1.1. Heuristics. Before giving more precise results below, we start with 
a heuristic discussion of why the relaxation time of a single card should 
fluctuate between order n 2 and order n?/ 2 (when k is of order n) according 
to the Diophantine approximations of k/n. We mention that Jonasson [7] 
already understood that the relaxation time is of order n 2 when k/n is near 
a simple rational, but conjectured that the relaxation time is of order nlogn 
for typical k. 
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Fig. 2. Shown here is the inverse spectral gap for the position of the <|kA for the n,k 
shuffle; the left picture plots the inverse gap as a function of k/n when n = 200 and k is 
varied from 1 to 199, while the right picture is for n = 1000 with k varying between 1 and 
999. There are spikes when k/n is well approximated by a simple rational number, and 
these spikes become more pronounced for larger n, with more of them visible. The predicted 
locations and heights of the peaks are shown in the n = 1000 picture. At the tops of the 
spikes the inverse spectral gap is 0(n 2 ), except for the spike at j where it is 0(n 3 ). At 
the low-lying areas between the spikes, the inverse spectral gap is 0(n 3 ^ 2 ). 



Suppose that the l|kA was initially at the top of the deck. Where will it 
be after t steps? Denote by s the number of visits of the l|kA to location 
1 (after time 0), by s the number of times it jumped from location n — k 
to 1, and by I its final location, which we assume to be less than n — k 
for simplicity. Clearly, s = 0(t/n) and s = s/2 ± 0(\fs). We get another 
approximate identity by noting that the card progresses at exactly unit 
speed from 1 to n — k and then either jumps to 1 or moves at approximately 
speed 1/2 from n — k to n. There are random fluctuations in how long the 
traversal from position n — k to n takes. We then have 

t = £- 1 + (n - k)s + 2k(s - 3) ± Q(^2k(s-s)). 

But it is likely that s — s~ t/(2n), and we assumed k is of order n, so this 
error term is Thus, we may rewrite this in the form 

£ = k(2s- s) - sn + t±e{Vi). 

In other words, at time t, the location £ G {1, . . . , n} of the l|kA is in one of 

2s - s = ±9( v / s) = ±®{\jt/n) 

intervals Ij of length @(y/t) each; these intervals are centered at the first 
@(y/t/n) multiples ±jk of k modulo n. A necessary condition for mixing is 
that these intervals cover most of {1, . . . , n}, which requires n = @(t/y/n), 
that is, t must be at least of order ra 3 / 2 . If the multiples {jk : \j\ < @(y/t/n)} 
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of k are close to uniformly distributed modulo n, and t > Cr?l 2 , then the 
intervals {Ij < G(y / t/n)} cover most of {1, . . . , n}. This happens when, 

for example, n/k is close to the golden ratio 1+ g , so in this case one might 
expect the relaxation time to be of order n 3 / 2 . On the other hand, if k/n is 
very close to a simple rational p/q, then there are only q distinct intervals 
Ij of length Q(y/i) each; for these intervals to cover most of {1, . . . , n}, the 
time t must be at least of order (n/q) 2 . 

1.2. Statement of results. We define gap(n, k) to be the spectral gap of 
a single card when the deck is shuffled according to the overlapping-cycles 
shuffle with parameters n and k. 



Theorem 1. Suppose thatp and q are integers such that < p/q < 1 and 
i for fixed p and q, when n is large enough and k = \ -° 

the spectral gap is 

gap(re,fc) = (l + o(l))- 



gcd(p, q) = 1. Then for fixed p and q, when n is large enough and k = [|n] , 

Tr 2 pq ( 1, ifp = qmod2, 



2n 2 \4, ifp^qmod2. 

It is interesting to consider the spectral gap when n is large but then held 
fixed while k is varied, as in Figure 2. We prove the following bounds: 



Theorem 2. 

(i) For any k G [1, re), we have gap(n, k) < (2ir 2 + o(l))n -3 / 2 . 

(ii) For any k G [l,re) , we have gap(n, k) > (4-7T 2 + o(l))n~ 3 . 

(hi) For almost all values of k G [1, n), u>e Ziaue £/iai gap (re, k) is of order 
n~ 3 / 2 . More precisely, when n is large enough, for any 5 > 0, for all but 
5 2 ^n values of k G [l,re), we have gap(?T,,fc) > (7r 2 /16 + o(l))5n~ 3 / 2 . 

Here the o(l) terms tend to as re — > oo independent of k or 5. 



The lower bound is tight, but as discussed in the open problem section, 
the upper bound can almost certainly be improved by a factor of y/2>. 

As the reader may have noticed already, the spikes in Figure 2 are actually 
shaped like narrow bells. These bells are shown more explicitly in Figure 3. 
Our next theorem characterizes the shape and width of these bells, by giving 
the spectral gap when k is close to but not exactly [(p/q)n \ . 

Theorem 3. Suppose that p and q are integers such that <p/q < 1 
and gcd(p, q) = 1. Then for fixed p and q, when n is large enough and 

\k — (p/q)n\ < re 3 / 4 x (constant depending on p and q), 
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Fig. 3. These plots show the same data from Figure 2 for the inverse spectral gap when 
n — 1000. On the left a specific bell curve is plotted on top of the spike at k/n = 1/3, and 
on the right such bell curves are plotted for all the simple fractions, including 2 and j . 



then the spectral gap is 

gap(n, k) = (1 + (D)^? (l + l l K l q) ) 

1, if p = q mod 2, 
4, if P^ q mod 2. 

When p/q = 0, the spectral gap is instead (1 + o(l))2ir 2 (k + k 2 )/n 3 , provided 
k<(n/2) 2 / 3 . 



Therefore, for each rational p/q ^ 0, the spectral gap ranges from order 
1/n 2 to order 1/n 3 / 2 when k/n is in a n _1//4 -neighborhood of p/q. 

For any given k and n, the spectral gap may be determined by the rational 
k/n or perhaps by a nearby rational with smaller numerator and denomina- 
tor. This leaves the question, for a given k and n, assuming n is large, how 
can we estimate the spectral gap? As Figure 3 indicates, the inverse spectral 
gap is determined by the rational whose bell is largest at k/n. 

Let us define 



(1) 



7r 2 / m 2 k 

j(n, k, m) = h 4n 

2n z V n 



m{n — k) 



2n 



7(n,/c)= min 7(77., k, m), 



where ||x|| is the number-theoretic norm ||x|| = dist(x,Z). 

Our main theorem on the spectral gap, from which the previous theorems 
can be derived as corollaries, is the following: 
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Theorem 4. For all large enough n, for any k such that 1 < k < n— 1, 
the spectral gap is 

gap(n, k) = (1 + o(l))j(n, k), 
where the o(l) goes to zero uniformly for the full range ofk, when n is large. 

This is a curious expression: there are terms in the min that do not cor- 
respond to eigenvalues, and there are eigenvalues that are not associated 
with terms in the min, but the sufficiently small terms do correspond to the 
eigenvalues with modulus sufficiently close to 1. Because the convergence 
holds uniformly for the full range of k, we can derive the previous theorems 
simply by studying 7(n,/c). For example, we have the following: 

Proof of Theorem 1. The min is achieved at m = q when p = q mod 
2, and m = 2q when p ^ q mod 2. □ 

An interesting consequence of Theorem 4 is the following: 

Theorem 5. For a G (0, 1), we have 

liminf7(n, |cm])ra 3 / 2 < 2ir 2 y/a/Vb, 
n — >oo 

and this bound is tight for a dense set o/a's. 

PROOF. Any irrational number x has infinitely many rational approx- 
imations p/q for which \x — p/q\ < l/(y / 5g |2 ), that is, \\qx\\ < l/(V5q) (see 
Hardy and Wright [5], Section 11.8). Even for rational x, there are still in- 
finitely many integers q for which \\qx\\ < l/{\fhq). Let us set x = (1 — a)/ 2 
and consider such an integer q for which q 2 > 1/a. We have, for n> q 2 , 

7 (n, [an] , q) = + 0(l/n)) + 4n||g[(l - a)/2\ + 0(q/n)\\ 2 ) 

it 2 (q 2 a 4nV 2 0(1] 



2n 3/2^ n l/2 5q 2 n l/2 

The optimal choice of n is n = \5q 4 a/4\ > q 2 , for which we get 
7 (n, L«nl,,)n 3 / 2 <^(^+ (l)), 

where the o(l) tends to as q — > oo, whence 

v -ci I "1\ 3/2 ^ 2-K 2 y/a 
hmrni7(n, \an \)n ' < 1= — . 
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There is a dense set of a's with continued fraction expansion that is even- 
tually all l's — these a's can be obtained from the golden ratio <f> by the 
Mobius transformation (acj) + b)/ (ccj) + d) , where a,b,c,d£Z and ad — be = 1 
([5], Section 10.11). For these a's, for any 5 > 0, there are only finitely 
many rational approximations p/q for which \\qx\\ < (1 — S)/(y/bq) ([5], Sec- 
tion 11.8), and so in the above calculations "<" may be replaced with 
= (1 + o(l))x, and so the above bound is tight. □ 

As the proof of the preceding theorem indicates, for a given k and n, 
the optimal choice of m for 7(n, k,m) will be the denominator of a rational 
approximation p/q to (1 — k/n)/2. \ipjq is not a continued fraction conver- 
gent, then there is a continued fraction convergent Pj/qj which both more 
closely approximates (1 — k/n)/2 and has smaller denominator qj < q, so 
j(n,k,qj) < 7(n, k,q). The best qj will never be larger than 0(n 1,/4 ), but it 
could be much smaller ifk/n has a particularly good rational approximation. 

Remark. Diaconis and Saloff-Coste [2, 3] analyzed the mixing times of 
two related shuffles. Goel [4] and Jonasson [7] analyzed a variant in which 
one of the bottom k cards is placed at the top of the deck. 

Roadmap. In Section 2 we review the characteristic polynomial for this 
Markov chain and make some preliminary observations about the approx- 
imate locations of the eigenvalues. There are many eigenvalues that are 
almost on the unit circle. Then in Section 3 we prove Theorem 4 charac- 
terizing the spectral gap. Using Theorem 4, we calculate in Section 4 the 
spectral gap for extremal and for typical values of k, proving Theorem 2. In 
Section 5 we deduce Theorem 3 for the case when k/n is well approximated 
by a simple rational. 

2. Eigenvalues for a single card. The Markov chain for a single card is 
as follows: A card in position i G [1, n — k — 1] moves to position i + 1. The 
card in position n — k moves to one of {l,n — k + 1}, each with probability 
i. A card in position i G {n — k + l,n — 1} moves to one of {i,i + 1}, and 
the card in position n moves to one of {l,n}, again each with probability 
^. Since most of the transition matrix is zero, it is easy to verify that its 
characteristic polynomial is 

\ n ~ k {\ - l/2) k - (A - l/2) k /2 - 2~ k ' 1 , 

and therefore, its eigenvalues satisfy 

(2) 2A n ^-l = (2A-l)~ fc . 

This equation for the eigenvalues first appeared in Jonasson's analysis [7], 
Section 4.1.2. 
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Fig. 4. Eigenvalues of the 50-20 shuffle. In general, the eigenvalues will be close to the 
unit circle near the (n — k)th roots of unity, or near the circle of radius 1/2 centered at 
1/2. Only eigenvalues close to 1 will be much closer to the unit circle than distance l/n. 



Proposition 6. The spectral gap for the position of a single card is 
positive. 

Proof. Note that if \z\ > 1, then \2z — 1| > 1, with equality iff z = 1. 
Thus, if |A| > 1, then the LHS of (2) is outside the unit circle and the RHS 
inside, and equality is possible iff A = 1. Since the derivative of (2X n ~ k — 
l)(2A-l) fc -l at A = l is 2n, the A = 1 eigenvalue has multiplicity 1. □ 

There are actually two reasons that we call the shuffle the overlapping- 
cycles shuffle. The first is that the cards are shuffled according to one of the 
two cyclic permutations (n,n — 1, . . . , 1) and (n — k, n — k — 1, . . . , 1). The 
other reason is that the eigenvalues for the position of a single card form two 
overlapping cycles, as can be seen from Figure 4. To see this, note that if A 
is sufficiently far inside the unit circle, then X n ~ k ~ and the LHS of (2) is 
approximately —1, so that 

(-1)1/* + 1 

where (— is any kth. root of —1 (sufficiently far from 1). If |A| ~ 1 but 
A is sufficiently far from 1, then (2A — l)~ fc « 0, so that (2) gives 

A»(l/2) 1/M) (l) 1/M) , 

where (l) 1 /( n_fc ) is any (n — /c)th root of 1 (sufficiently far from 1). The most 
interesting range is where A ~ 1, that is, where the two cycles overlap, and 
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as we shall see in the next section, it is these eigenvalues that will ultimately 
determine the spectral gap of the Markov chain for a single card. 

3. Spectral gap for a single card. We proceed to analyze the eigenvalues 
using equation (2) when n is large. For an eigenvalue A of the overlapping 
cycles shuffle with parameters n and k, we shall define e, a and b to satisfy 

A = exp[— e + ma\ where— 1 < a < 1, and 

A n B = exp[— (n — k)£ + inb] where — 1 < 6 < 1. 

At this point we restrict our attention to eigenvalues for which e = o(l/n). 
We shall see that there are always eigenvalues for which e = 0(ra~ 3 / 2 ), and 
for some values of n and k there are eigenvalues with much smaller e. We 
start with a lemma showing that such eigenvalues satisfy A ~ 1 and X n ~ k ~ 1, 
and derive successively more precise estimates of A. 

Lemma 7. With the notation (3), for any eigenvalue A for which e = 
o(l/n), we have a = o(l) and b = o(l). 

Proof. We have |2A — 1| > 1 — 2e + c(a), where c(a) > for a / is a 
positive constant depending on a. Then 

|2A n ^ fc | - 1 < |2A n " fc - 1| = |2A - l\~ k < (1 - 2e + c(o))" fc . 

When e = o(l) and a = 6(1), it follows that |A n_fc | < 1 — c(a), where c(a) is 
a positive constant depending on a, which contradicts e = o(l/n). 

We also have \2\ n ~ k — 1| > 1 — ne + c{b) where c(b) > is as above. Thus, 

|2A| - 1 < |2A - 1| = |2A n ~ fc - < (1 - n£ + c{b)y 1/k . 

When £ = o(l/n) and b = 6(1), it follows that |A| < 1 — c(b)/k, where c(b) 
is a positive constant depending on b, which contradicts e = o(l/n). □ 

Lemma 8. With the notation (3), for any eigenvalue A for which e = 
o(l/n), 

|ott| < ^2n£/k + 0{{n£/k) 3/2 ), 
H <y2^i + 0((?ie) 3/2 ). 

Proof. Let z = e _:r+ * 6 ' where a; and # are real and x > 0. By the triangle 
inequality, \ + \z — ^| > \z\ > 1 — x, so 

(4) |2z-l|>l-2x. 
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But by the law of cosines, 

\2z- 1| 2 = l + A\z\ 2 - 4|^| cos6» 

(5) = 1 + 4 - 8x + 0{x 2 ) - 4(1 - x + 0(x 2 )){l - 9 2 /2 + 0(6> 4 )) 
= l-Ax + 29 2 + 0(x 2 ) + 0(x9 2 ) + 0{9 4 ). 

From (2), using (4) with z = A, we get 

(6) |2A n " fc - 1| 2 = |(2A - iy k \ 2 < (1 - 2e)~ 2k = 1 + 4fe + 0(n 2 e 2 ), 
and using (5) with z = X n ~ k , we get 

(7) |2A n " fc - 1| 2 = 1 - 4(n - fc)e + 0(n 2 e 2 ) + 2&V(1 + O(ne) + 0(b 2 )), 
and upon combining (6) with (7), we get 

4rae + 0(n 2 e 2 ) > 2&V(1 + 0(ne) + 0(6 2 )). 

Hence, taking square roots, 

V2ns x (l + 0(ne)) > |6tt|(1 + O(ne) + 0(6 2 )). 

Since e = o(l/n), and 6 2 = o(l) by Lemma 7, we learn b = 0{^/ne), and 
hence, 

|6vr| <V2^ex (l + 0(ne)). 
From (2), using (4) with z = X n ~ k , we get 

|2A- l| 2 = |2A n - fc -ir 2/fe 

< (l-2(n-k)sy 2/k 

(8) 

= i + ^ + o(KAf), 

and using (5) with z = A, we get 

(9) |2A - 1| 2 = 1 - 4e + 0(e 2 ) + 2aV(l + 0(e) + 0(a 2 )), 

and upon combining (8) with (9), we get 

Ane/k + 0({ne/k) 2 ) + 0{e 2 ) > 2aV(l + 0(e) + 0(a 2 )). 
Hence, taking square roots, 

^2ne/k x (1 + 0(ne/k)) > |avr|(l + 0(e) + 0(a 2 )). 

Since e = o(l/n), and a 2 = o(l) by Lemma 7, we learn |a| = 0(y/ne/k), and 
hence, 

M < ^2ne/k x (1 + 0(ne/k)). □ 
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Lemma 9. With the notation (3), for any eigenvalue A for which e = 
o(l/n), we can write 

IT 2 



nim/n (k(m/n) + (m(l — k/n) mod 2) ) 

2n 



X = exp 

(10) 



for some integer m = (9(n 3//2 e 1 / 2 /A; 1 / 2 ) = o{n/k 1 / 2 ) , where mod 2 is inter- 
preted so that — 1 < * mod 2 < 1. 

Proof. We apply the following expansion for z ~ 0: 
(11) 2e z -l = exp[2z-z 2 + 0(z 3 )]. 

By Lemma 8, we know that eigenvalues A for which e = o(l/n) satisfy 
| logA n_fc | = 0(y/ne), so 

2\ n ~ k - 1 = exp[-2(n - k)e + 2biti + + 0((ne) 3 / 2 )]. 



Lemma 8 also tells us | log A| < 0{^en/k), so 

2A - 1 = exp[-2e + 2am + a 2 ir 2 + 0{{en/kf/ 2 )}, 
(2A - l)" fc = exp[2A;e - ka2m - ka 2 ir 2 + 0((en) 3/2 /k 1/2 )}. 

For any eigenvalue A, by (2), these expansions are equal: 
-2(n - k)e + 2bni + b 2 n 2 

(12) 

= 2ks - ka2iri - ka 2 iT 2 + 0((ne) 3/2 ) mod 2ni. 
Comparing the real parts of (12), we obtain 

- 2ne + &V = -ka 2 ir 2 + 0((ne) 3/2 ) 

(13) e = l!( / fca 2 + 6 2 ) + 0(n 1 /V/ 2 ). 

2n 

Comparing the imaginary parts of (12), we obtain 

b = -ka + 0{{nef /2 ) mod 1. 
Of course, from (3) it is obvious that 

b = (n — k)a mod 2. 
Subtracting these two relations, we obtain 

= na + O{(nef /2 ) mod 1. 
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Since na is very close to an integer, we define 

m = [na] = 0(n 3/2 e 1/2 /k 1 ' 2 ) = o{n/k 1 ' 2 ). 

Then 

a = m/n + 0(e 3/2 n 1/2 ), 
b = [m(l - k/n) + 0(e 3/2 n 3/2 )] mod z, 
where mod 2 is interpreted so that — 1 < * mod 2 < 1. We have then that 
A = exp[a7ri — e] 

7T 2 

= exp irim/n (k(m/n) 2 + (mil — k/n) mod 2) 2 ) + 0(n l ^ 2 e 3 ^ 2 ) 

[ In 

This estimate is valid for any eigenvalue A for which e = o(l/n). □ 

Lemma 10. If m is an integer such that \m\ =o{n/^/k) and |m(l — 
k/n) mod 2| = oil), then there is an eigenvalue of the form given in equa- 
tion (10). 

PROOF. Let us define 

b' = m(l - k/n) mod 2 with -1 < b' < 1, 

X = ^(k(m/n) 2 + b' 2 ), 
in 

X(z) = exp[irim/n — X + z], 

f(z) = (2X(z)-l) k (2X(z) n - k -l)-l. 

Our aim is to show that there is some choice of z S C such that f{z) = 
and \z\ = Oin l / 2 X^/ 2 ). Of course, if f(z) = 0, then A (z) is an eigenvalue. 
Our hypotheses imply X = o(l/n), so \z\ = o(X), e = (1 + o(l))X, and \z\ = 
0(n 1 / 2 e 3 / 2 ), so the existence of such a z will prove the lemma. 
Since 

\(z) n - k = exp[7ri&' -(n-k)(X- z)}, 

the expansion (11) together with \irb'\ < \JlnX = o(l) imply that whenever 
\z\ = 0(X) we have 

2X(z) n - k - 1 = exp[27ri6' - 2(n - k){X - z) + ir 2 b' 2 + 0(n 3 / 2 X 3 / 2 )]. 



The expansion (11) together with \irm/n\ < \J1nX/k = oil) imply 

2\(z) - 1 = exp[2vrim/n - 2{X - z) + vr 2 (m/n) 2 + 0((Xn/k) 3/2 )\. 
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Thus, modulo 2iri, 

log(l + f(z)) = k[2mm/n - 2(X - z) + ir 2 (m/n) 2 + 0((Xn/k) 3/2 )] 
+ [2mb' - 2(n - k)(X - z) + vr 2 6 /2 + 0{n 3 / 2 X 3 / 2 )) 
= 2irikm/n + ir 2 k(m/n) 2 + 2-Kib 1 - 2n(X - z) + ir 2 b' 2 

+ 0(n 3 ' 2 X 3 l 2 ) 
= 2nz + 2-Kikm/n + 2mb' + 0{n 3 / 2 X 3 / 2 ) 
= 2nz + 0(n 3/2 X 3/2 ) (mod 2m), 

whence 

f(z) = 2nz + 0{n 3/2 X 3/2 ). 

By considering z's for which \z\ is a large enough constant times n l / 2 X 3 / 2 
[but still within the constraint \z\ = O(X)], we may apply Rouche's theorem 
(see [1], Section 4.5.2) to deduce that there is a z for which \z\=0{n 1 / 2 X 3 / 2 ) 
and f (z) =0. □ 

Recall the definition of 7(71, k) from equation (1). We are now ready to 
prove Theorem 4, that is, that the spectral gap is gap(n, k) = (l + o(l))j(n, k). 
We actually show a bit more, that gap(n, k) = j(n, k) + 0(n 1 / 2 7(n, k) 3 / 2 ). 

PROOF of Theorem 4. Recall the definition 

7T 2 

7(71, fc, m ) := — (kim/n) 2 + (m(l — kin) mod 2) ), 
2n 

which is the subexpression corresponding to e (ignoring the error terms) in 
equation (10). 

For any real numbers a and N > 1, there is an integer m so that 1 < m < 
iV and Track is within l/([A r J + 1) < 1/N of an integer (see [5], Section 3.8, 
Theorem 36). For our purposes, we are interested in a = ^(1 — k/n) and 
N = 2 1 / 2 ?7 1 / 2 //c 1//4 . Then there is a positive integer m < N for which 

7 (ra, fe,m) < — (k(N/n) 2 + (2/iV) 2 ) = ^— < 2vr 2 ra- 3 / 2 . 

2ra n z 

Let 7ra denote the positive integer minimizing 7(72, k, rh), that is, 7(77, k, rh) = 
7(77, k). We have 

— k(m n) <7(ra,K,7ra)< ^ — , 

2ra 

so |rn| < 2n l / 2 k~ 1 ^ = o(n/k 1 / 2 ), and similarly, 

— (771(1 — k/n) mod 2) 2 < 7(77, k, m) < ^ — > 
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so \fh(l — k/n) mod 2| < 2k 1 ^n x l 2 = o(l). Since to satisfies the hypotheses 
of Lemma 10, it gives us an eigenvalue A for which e = (1 + o(l))7(n, k, fa) = 
(l + o(l))7(n,A:). 

Since 7(71, fc) < 2-7r 2 n 3//2 , this proves Theorem 2 part (i). 

Therefore, any eigenvalue A which achieves the spectral gap will satisfy 
e = o(l/n), so we can use Lemma 9 to express it according to equation (10) 
for some integer to = to*. The spectral gap is then (1 + o(l))7(n, k, to*). 
We claim further that to* 7^ 0, since otherwise A = exp so e = 

0(n l/2 £ 3/2 )=o(e)) 

so e = and A = 1 , which by Proposition 6 is not the 
eigenvalue achieving the spectral gap. Since there is an eigenvalue for which 
e = (1 + o(l))7(re, k) < (1 + o(l))7(n, A;, to*), we conclude that the spectral 
gap is in fact (1 + o(l))7(n, k). □ 

4. Spectral gap for extremal and typical values of k. Recall that the 
upper bound on the spectral gap [Theorem 2 part (i)] was shown in the 
proof of Theorem 4. We proceed to prove the lower bound. 

Proof of Theorem 2 part (ii). If for some nonzero k and m we have 
7(77., k, to) <4-7r 2 n -3 , then 

7T 2 

— k(m/n) 2 < 7(71, k, m) < Air 2 n 3 , 
2n 

and so km 2 < 8. Similarly, 



r 2 

-(m(l — k/n) mod 2) < 7(74, k, to) < 4-7r 2 n~ 



2n 

so |to(1 — k/n) mod 2| < y/8/n. Since k < 8, for n > 11 this implies that 
1 771 1 7^ 1, and hence, m = ±2 and k < 2. Then for n > 11, 

7T 2 27T 2 

7 (n, fe,m) = — (k{2/n) 2 + (2A;/n) 2 ) = — T (k + fc 2 ), 

and our assumption implies k = 1 and j(n,k,m) = 47r 2 n -3 . Thus, whenever 
ri > 11 we have 7(77, k) > 47r 2 n~ 3 , and so the result follows from Theorem 4. 
□ 

We now prove the last part of Theorem 2 which pertains to "typical" 
values of k. 

Proof of Theorem 2 part (iii). We will show that for any n > 1 and 
any 6 > 0, for all but 4^ 2 / 3 n values of k, we have 7(77, k) > (tt 2 /2)8n~ 3 ^ 2 . 
Let us start by assuming k > b 2 ^n. If to, > <5 1 /6 n i/4 ) then 

. Tr 2 m 2 k tt 2 5 

7(77, fc, to) > - > 



2n 3 - 2n 3 / 2 ' 
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We also have 

7T 2 7T 2 

7(77,, fc,m) > — (m(l — k/n) mod 2) 2 = — ^(mk* mod 2n) 2 , 
2n 2n 6 

where k* =n — k and mod denotes the centered mod. If k* is a uniformly 
random integer from the range 1 < k* < 2n, then mk* mod 2n is uniformly 
distributed among the multiples of gcd(m,2n). For any x with < x < n, 
there are 



1 gcd(m, 2n) < 2x + gcd(?n, 2n) 
gcd(m, 2n) J 

values of k* (1 < k* < 2n) for which \mk* mod 2n| < x, and a fortiori at most 
x + ^ gcd(?n, 2n) values of A;* (1 < k* < n — 1) for which \mk* mod 2n| < x. 
Upon taking x = <5 1//2 n 3 / 4 , and summing over 1 < m < <5 1//6 n 1//4 , we find that 
there are at most 

y j I x + -gcd(m,2?7,J J < o ' n-\ 

m=l 

values of k > 5 2 ^ 3 n and m < 5 1 / 6 n 1//4 for which 7(n,/c,m) < (7r 2 /2)<5n~ 3 / 2 . 
Thus, there are at most 

(14) 25 2 / 3 n + ^ 1 / 3 n 1 / 2 + ^ 1 / 6 n 1 / 4 

values of k for which 7(72, /c) < (7r 2 /2)<5n -3 / 2 . 

If <5 2//3 n < I, then the bound (14) is less than one, so there are < 
4£ 2 / 3 n values of k for which j(n,k) < {it 2 /2)<5n~ 3 / 2 . If 5 2 / 3 n > \, then the 
bound (14) is < 45 2 / 3 n. □ 

5. The narrow bells near simple rational k/n. Let us suppose that 
k/n is well approximated by a simple rational p/q where <p/q< 1 and 
gcd(p, q) = 1. Our aim here is to prove Theorem 3, which assumes that p 
and q are held fixed while n — > 00. 

Proof of Theorem 3. Let c be defined by k = (p/q)n + cn 3 / 4 . Recall 

2 

7(ra, fe, m) = — ^(m 2 k/n + n(m(l — k/n) mod 2) 2 ). 
2n^ 

Consider the term 

n(m(l — fc/n) mod 2) 2 = n(m(l — p/g — cn~ 1 ^) mod 2) 2 . 

As we saw in the proof of Theorem 1, if c = 0, this term is minimized at 
m = Aq, where 

1, if p = q mod 2, 



,4: 



2, if p ^ g mod 2. 
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When m is divisible by Aq and men 1 / 4 < 1 we have 

7T 2 

(15) 7(n, k, m) = -^—^m 2 {k/n + c 2 n x l 2 ). 

If m is not divisible by Aq, then men -1 / 4 would have to be as large as 
(1 — o(l))2/(Aq) before the second term of j(n, k, m) would be smaller than 
0(l/n). But for m this large, we have 

, , . 7T 2 k An 1 / 2 

(16) J{n,k,m)>— — 



2n 2 n A 2 q 2 c 2 

on account of the first term. 

If c is small enough for the first term of (15) to be nonnegligible, then 
(16) is order 1/n and thus irrelevant to computation of 7 (n, k). If the first 
term of (15) is negligible, then we may compare (15) to (16), and we find 
that (15) determines 7 (n,A;), provided 

A 2 q 2 c 2 < k 4 

~~ n A 2 q 2 c 2 ' 

< (Ak/n) 1 / 4 
C " Aq ' 

If p + and c < {Ap/q) l / 4 /{Aq), we obtain 

7 (n, k) = (1 + o(l))^A 2 q 2 (p/q + cV/ 2 ). 
When p = we instead obtain 

2vr 2 
n 

for k< (n/2) 2 / 3 . □ 

6. Open problems. 

1. We would like to understand the mixing time of the over lapping-cycles 
shuffle as a Markov chain on the symmetric group. Is it within a factor 
of O(logra) of the relaxation time of a single card? 

2. For all large n there is a A; of the form k = n — (2/3 1 / 4 + o(l))n 1 / 4 
such that 7(71, k) = (2tt 2 /\^3 + o(l))n _3//2 . In fact, it appears that for 
each real a £ (0,1] there is a /c of the form k = (a + o(l))n for which 
7(n, A;) = (2tt 2 /\/3 + o(l))fe 1 / 2 n -2 . Furthermore, it appears that 7(71, A;) 
can be no larger than this, as is illustrated in Figure 5. Does the inequality 
gap(n,A;) < (2tt 2 />/3 + o(l))fc 1 /2 n -2 hold? 



7 (n,A;) = (l + (l)) — (fc + A: 2 ) 
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Fig. 5. Shown here is the relaxation time for the position of the 4k A for the overlap- 
ping-cycles shuffle with parameters n = 1000 and k, together with a lower envelope on this 
function coming from a hypothetical bound "f(n, k) < (2-7r 2 /\/3 + o(l))fc 1 ^ 2 n -2 . 
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